Entropic transport of finite size particles 
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Transport of spherical Brownian particles of finite size possessing radii R < i? max 
through narrow channels with varying cross-section area is considered. Applying the 
so-called Fick-Jacobs approximation, i.e. assuming fast equilibration in orthogonal 
direction of the channel axis, the 2D problem can be described by a ID effective 
dynamics in which bottlenecks cause entropic barriers. Geometrical confinements 
result in entropic barriers which the particles have to overcome in order to proceed in 
transport direction. The analytic findings for the nonlinear mobility for the transport 
are compared with precise numerical simulation results. The dependence of the 
nonlinear mobility on the particle size exhibits a striking resonance-like behavior 
as a function of the relative particle size p = R/R m3iX ; this latter feature renders 
possible new effective particle separation scenarios. 

I. INTRODUCTION 

The diffusive behavior of Brownian particles depends mainly on their size, the interaction 
between them, and the environment where they are situated in. If, in addition to these 
characteristics, particles are confined within narrow, tortuous structures such as nanopores, 
zeolites, biological cells and microfluidic devices, the restriction of the space available for the 
particles will cause entropic barriers that will have strong impact on the diffusive behavior 
(cf. Ref. []J and references therein). Effective control schemes for transport in these systems 
require a detailed understanding of the diffusive mechanisms involving small objects and, 
in this regard, an operative measure to gauge the role of fluctuations. The study of these 
transport phenomena is, in many respects, equivalent to an investigation of geometrical 
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constrained Brownian dynamics. As the role of inertia for the motion of the particles through 
these structures can typically be neglected, the Brownian dynamics can safely be analyzed 
by solving the Smoluchowski equation in the domain defined by the available free space upon 
imposing reflecting boundary conditions at the domain walls. 

However, solving the boundary problem in the case of nontrivial, corrugated domains 
presents a difficult task. A way to circumvent this difficulty consists in coarsening the de- 
scription by reducing the dimensionality of the system considering only the main transport 
direction, but taking into account the physically available space by means of an entropic 
potential [l 3 ] . The resulting kinetic equation for the probability distribution, the so-called 
Fick- Jacobs equation, is similar in form to the Smoluchowski equation, but contains now en- 
tropic contributions leading to genuine dynamics which distinctly differs from those observed 
for purely energetic potentials. 

The driven transport of particles across bottlenecks [l-4], such as ion transport through 
artificial nanopores or artificial ion pumps js-Sj or in biological channels js~ 12] , are striking 
examples where the diffusive transport is regulated by entropic barriers. In addition, geo- 



metrical confinements and entropic 
Stochastic Resonance phenomenon 



3arriers play also a prominent role in the context of the 



13|4l7|. 



Our objective with this work is to investigate the mobility of noninteracting spherical 
Brownian particles in channels with varying cross-section width. In particular, we are inter- 
ested in the influence of the particle size on the transport within a periodic entropic potential 
exhibiting barriers which arise from the geometrical restrictions. 

The paper is organized as follows: in section 2 we introduce the model and define the 
theoretical and numerical problem. Further on, in section 3 we present the basic principles 
of the Fick- Jacobs approximation allowing for reducing the two-dimensional problem to an 
one-dimensional one. The results are presented in section 4. Finally, we give the main 
conclusions in section 5. 



II. MODELLING 

Transport through pores or channels (like the one depicted in Fig. [1]) may be caused by 
different particle concentrations maintained at the ends of the channel, or by the application 
of external forces acting on the particles. Here, we exclusively consider the case of force 
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driven transport of spherical particles of radius R. The external force F = Fe x is pointing 
parallel to the direction of the channel axis. As small deviation from this assumption does not 
affect our results, certainly not within the limits of validity of the Ficks- Jacob approximation. 
Moreover, we shall assume low concentrations of spherical particles such that particle-particle 
interactions and all hydrodynamic interaction effects can consistently be neglected. 



A. Dynamics inside the channel 

In general the dynamics of a suspended Brownian particle is overdamped [18( and well 
described by the Langevin equation: 

r lR% = F+ y/vRksmt), (1) 

where f denotes the center position of the spherical particle in the two-dimensional channel, 
k-Q the Boltzmann constant, T the temperature and is the standard 2D Gaussian noise 
with (£(£)) = and (£i (£)£/(£')) = 25y-5(£ — t') for i,j = x,y. The friction coefficient tjr is 
given by Stokes' law: 

Vr = 6nvR (2) 

and depends on the shear viscosity v of the fluid and the particle radius R. In addition to 
Eq. [1] the full problem is set up by imposing reflecting boundary conditions at the channel 
walls. The boundary of the 2D periodic channel which is mirror symmetric about its x-axis 
is given by the periodic function y = ±u(x) with u(x + L) = u(x) where L is the periodicity 
of the channel. oo m i n and w max refer to the half of the maximum and minimum channel width, 
respectively. 

To further simplify the treatment of this problem, we introduce dimensionless variables. 
We measure all lengths in units of the periodicity of the channel, i.e. x = x'L. As a unit of 
time t we choose twice the time it takes the largest transporting particle to diffusively cover 
the distance L which is given by r = L 2 r] max /(k-BT), hence t = rt' . The largest transportable 
particle is the particle with radius i? max = ^Vm- Accordingly, the friction coefficient of a 
particle of radius R is then given by rj = pi] max with the ratio of spherical particle radii 
being p = Rj R mSkX and r/ max = 6nuR max . 

Summarizing, the Langevin equation ([T|) reads in dimensionless variables: 

df* f fl - 
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FIG. 1: (Color online) Sketch of the 2D periodic channel with periodicity L, the minimum half 
channel width uj m i n and the maximum half channel width u; max . The spherical Brownian particle 
of radius R is subjected to the force F. 



where the dimensionless force parameter 
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For the sake of better readability, we shall skip all the primes in the following and proceed, 
if not mentioned explicitly otherwise, with dimensionless variables. 

The corresponding Fokker-Planck equation for the time evolution of the probability dis- 



tribution P(r,t) takes the form 19] 



Of 

where J(r,t) is the probability current: 

J(r,t) = - p (/e;-v)p(f,t). (6) 

B. Boundary conditions 

As the particles are confined by the channel structure, the probability current has to 
vanish at the boundaries. Due to the finite size of the particles their center position can 
aproach the boundary only up to its radius. Consequently, the position vector r of a particle 
with radius R never approaches the channel walls and is restricted to only a portion of 
the inner channel area, cf. Fig. [TJ The effective boundary function u eS (x), which serves as 
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FIG. 2: (Color online) Sketch of the original tube geometry given by the boundary function ui(x) 
and the effective boundary function u eS (x) for the center of the spherical particle. Hereby, the 
effective boundary function depends on the radius of the particle which is given in dimensionless 
units by p: p = 0.27 (red dotted line), p = 0.81 (blue dashed line). 

boundary for the center of mass, exhibits the distance R from the original, true boundary 
function u(x). Consequently, the "no-flow" boundary conditions for the center of mass 
dynamics read: 

J(f,t)-n = 0, for f 6 effective boundaries, (7) 

where n denotes the normal vector field at the effective channel walls. For the considered 
2D channel structure, the boundary condition becomes 

^{ /i>(liVi() _a^} + a^ = o, (8) 

at y — ±cu eS (x). 

Note, that the effective boundary function exhibits a complex dependence on the particle's 
radius R and could not be given explicitly. If the curvature of the channel wall function 
oj{x) is larger than that of the particle, the effective boundary function exhibits a kink, cf. 
Fig. 13 

For an arbitrary form of u)(x), the boundary value problem defined by Eqs. (jSJ), ([6]) and 
(jHJ) is very difficult to solve. Despite the inherent complexity of this problem an approxi- 
mate solution can be found by introducing an effective one-dimensional description where 



geometric constraints and bottlenecks are considered as entropic barriers j^, ^. I20M25 ] . 
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III. FICK-JACOBS APPROXIMATION 

The ID equation is obtained from the full 2D Smoluchowski equation upon the elimina- 
tion of the transversal y coordinate assuming fast equilibration in the transversal channel 
direction. 

A. The Fick-Jacobs equation 

The marginal probability density along the axis of the channel is defined by 

/uj cfi (x) 
P(x,y,t)dy. (9) 
■ui cS {x) 

Assuming fast equilibration in y-direction the 2D probability distribution becomes 

P(x,y,t) = V(x,t)Q(y\x), (10) 

with the local equilibrium distribution Q(y\x) of y, conditional on a given x. If there is 
no force component in transversal channel direction (as it is in our case), the conditional 
distribution Q(y\x) is uniform and reads due to the normalization condition: 

Q(y\x) = l/(2u cS (x)). (11) 

Then on integrating the full 2D Smoluchowski equation ([5]) and making use of Eqs. (Q, 

ffTUj) and fTTT]) . the Fick-Jacobs equation for the spherical particle is obtained: 

dV(x,t) _ 1 d f dA(x) d_\ v( f) (u) 
dx pdx \ dx dx) 

with the dimensionless free energy A(x) = —fx — lnu; (x). For a periodic channel this 

free energy assumes the form of a tilted periodic potential with the bottlenecks forming 

entropic potential barriers. Note, that for a straight channel, i.e. constant effective boundary 

function, the entropic contribution vanishes and the particle is solely driven by the external 

force. 

Introducing the x-dependent diffusion coefficient D(x) in Eq. [12] considerably improves 



the accuracy of the kinetic equation, extending its validity to more winding structures 
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251 ] . The expression for D(x) (in dimensionless units) 



D(x) = - (13) 

[1 + {du{x)/dxf] ' 



has been shown to appropriately account for curvature effects of the confining walls 
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B. Nonlinear Mobility 

Besides the effective diffusion coefficient, the average particle current, or equivalently the 
nonlinear mobility serves as key quantity of particle transport through periodic channels. 
For any non-negative force the average particle current in periodic structures can be obtained 



from Ref. 
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(x) = (t(x ^x + l))- 1 , (14) 

where (t(a — > b)) denotes the mean first passage time of particles starting at x = a to 
arrive at x = b. Within the Fick- Jacobs equation (fl2j) . the mean first passage time can be 
determined: 

(t(a -+b))=p f ' dxexp(-fx)/co cS (x) f dyexp(f y)uj cS (y) . (15) 

J a J — oo 

The nonlinear mobility p(f) is defined by p(f) = (x) / f and can be obtained as 

1 - cxp(-/) 

/ / dzl(z, f) 
Jo 



1 1 - exp(-/) 
M/) = - 7i > ( 16 ) 



where 

fX 



I(zJ) = exp(-fx)/u cS (x) [ dyex P (fy)cu cS (y) . (17) 

Jx-l 

In case of a straight channel with u m i n = cu max , an exact analytical solution of the full 
2D Smoluchowski equation ([5]) is known and the nonlinear mobility equals the free mobility 
(i.e. without geometrical constrictions) 

= y^frcc = 1/P = -Rmax/-R (for straight channels) . (18) 

Consequently, the influence of the confinement can be expressed by the ratio of nonlinear 
mobility for the transport through the channel and the one for the unrestricted case: 

M/) = 1 - exp(-/) ^ 



JO 

IV. PRECISE NUMERICS FOR A TWO-DIMENSIONAL CHANNEL 

GEOMETRY 



The nonlinear mobility, predicted analytically within the Fick- Jacobs approximation, has 
been compared with Brownian dynamic simulations performed by a numerical integration 
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of the full 2D Langevin equation ([3]), using the stochastic Euler algorithm. As random 
number generator we used the Box-Muller- and MT19937-algorithm from the GSL library. 
The sinusoidal shape of the considered two-dimensional channel is described by 



u(x) := a sin(27rx) + b , 



(20) 



with the two dimensionless channel parameters a and b. In physical units, these two pa- 
rameters are given by aL and bL, respectively. Note that u(x) may also be regarded as the 
first terms of the Fourier series of a more complex boundary function. Due to the symmetry 
with respect to the x-axis, the boundary function could be given in terms of the maximum 
half-width w max = b + a and the aspect ratio of minimum and maximum channel width 



C = ^min/uVax (with UJ n 



aj, i.e. 



UJ(X) 



uj(x) 



;i-e) 



sin (2t[x) + 



sin (27rx) 



+ w n 



1 + e 



(21) 
(22) 



To ensure, that the spherical particles of radius p stay within this channel geometry, the 
integration was carried out performing "no-flow" boundary conditions at the channel walls. 
By averaging over 10 5 simulations we obtain the steady-state average particle current 

(*(*)> 



(x) 



lim 

t— ¥00 



(23) 



and the nonlinear mobility p = (x)/f. 



A. Nonlinear mobility: force and temperature - dependence 



Figure [3] depicts the nonlinear mobility as a function of the scaling parameter / for two 
different particle radii and a fixed channel geometry: u(x) = 0.7/(2-71") sin(27rs) + 1.02/(2-71"). 
Strikingly, the transport through such channel structures is distinctly different from the one 



29j . This phenomenon is due 



occurring in one-dimensional periodic energetic potentials |2|-|4|, 
to the different temperature dependence of the barrier shapes. Decreasing the temperature in 
an energetic periodic potential decreases the transition rates from one cell to the neighboring 
one by decreasing the Arrhenius factor 30] and, therefore, reduces the nonlinear mobility. 
For the periodic channel system, a decrease of temperature results in an increase of the 
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FIG. 3: (Color online) Graph for the scaled nonlinear mobility as function of the force parameter /. 
In the Langevin simulation the different symbols correspond to different particle radius of p = 0.2 
(red triangles) and p = 0.8 (blue diamonds). The relative error of the simulation results is smaller 
than 0.01. The Fick-Jacobs results, Eq. (|19p . correspond to the solid lines. The boundary function 
reads: co(x) = (0.7/2vr) sin(27nr) + 1.02/2vr. 



dimensionless force parameter /, cf. Eq. (jlj) and consequently, in a monotonic increase of 
the nonlinear mobility. 

Due the geometrical restrictions, the nonlinear mobility is always smaller than the mo- 
bility for the free case, cf. Fig. [3j With increasing scaling parameter, the nonlinear mobility 
tends to that of the free case, i.e. \i — > /Xf ree for / — > oo. 

A comparison of the analytics obtained by means of the Fick- Jacobs-approximation with 
the precise numerics enables one to determine validity criteria for the Fick-Jacobs approx- 
imation, for further details see Ref. jj], 4]. According to them the applicability of the 
Fick-Jacobs approximation for the transport of point particles depends on the smoothness 
of the geometry and the scaling parameter /. For finite size particles, the diameter of the 
particles should become an additional parameter in the validity criteria. In particular, the 
particle's radius determines the maximum width of the effective channel structure. A larger 
particle leads to a smaller maximum effective width. Since for a fast equilibration in the 
orthogonal tube direction the timescale for the orthogonal diffusion process must be smaller 
than the timescale for the drift j^], a smaller maximum channel width favors the valid- 
ity of the Fick-Jacobs approximation. Thus, with increasing particle size the range of the 
applicability of FJ increases, cf. Fig. [3j 
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However, it turned out, that the transport phenomena presented below, occur for channel 
structures for which the validity criteria is not fulfilled. Therefore, we stick in the following, 
to the numerical results only. 

B. Particle size 

Surely, the transport of particles through small channel systems depends on the size of 
the particles. In particular, the effect of the size on the nonlinear mobility is two-fold. 
Firstly, the friction coefficient depends on the particle size, resulting in a p-dependence of 
the nonlinear mobility even for the case of unconstrained motion (free case), cf. Eq. ffT8|) . 
With increasing radii ratio p, the mobility declines. Secondly, there is the influence of the 
geometrical confinement. The extent of the bottleneck (2o; wm i n ) gives a limit to the size 
of particles able to travel through the channel. In our scaling the largest sphere possible 
to overcome the geometry's bottleneck has a radii ratio of p = 1. Considering particles of 
different sizes, the effective bottleneck (2u^ in ) will be smaller for spheres of higher diameter. 
It is intuitive that a small bottleneck hinders the transport. 

Figure Hl^a) depicts the nonlinear mobility as function of the radii ratio p for different 
geometries. For a straight channel one observes the nonlinear mobility of the free case /if ree 
which depends reciprocally on the radii ratio p. In presence of geometrical restrictions, i.e. 
for varying cross-section width, the nonlinear mobility is smaller than /if ree , cf. Fig. H] (a). 
Moreover, upon decreasing the bottleneck half-width (i.e. decreasing the aspect ratio e) of 
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the structure, the nonlinear mobility decreases 

Deviations from the l/p- dependence show another effect of the geometrical confinement, 
cf. Fig. H] (a). In order to focus on the geometrical effect, we consider the scaled nonlinear 
mobility, i.e. the nonlinear mobility relative to the nonlinear mobility in free case: p/p.f Tee , 
cf. Fig. H] (b). This is equivalent to the consideration of the nonlinear mobility of a point 
particle moving in the effective channel geometry defined by the effective boundary function 
cu (x), which still depends on the parameter p. With increasing p, the maximum half- width 
of the effective geometry shrinks. As a consequence, the sojourn time, the particle spends 
on average in a bulge of the channel structure decreases with increasing p and the mobility 
of the point particle in the effective geometry grows. This behavior causes the maximum in 
the scaled nonlinear mobility and the shoulder in the dependence of the nonlinear mobility 
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FIG. 4: (Color online) The numerically obtained non-linear mobility \x (a) and the scaled nonlinear 
mobility p/pt ree (b) as function of the radii ratio p = R/Rmax for different channel geometries 
(constant-width-scaling: cj max = const), cf. Eq. (f22j) . 

on the radius, cf. Fig. HI 

C. Role of the channel structure 



The confinement by the considered channel geometry can be altered by systematically 
changing the parameters Co> max and co m i Q or w max and e in the boundary function, Eq. (I2T]) or 
Eq. (1221) respectively. While for Fig. H] we examined a constant maximum half- width w max 
= const and varied the aspect ratio e which is equivalent to vary the half-width u m i Q at 
the bottleneck, cf. the constant-width-scaling in Ref. 3l|, it is instructive to also consider 
both, the constant-bottleneck-scaling u m i n = const as well as the constant-ratio-scaling e = 

W m in/W max = COnst. 

As we keep the bottleneck-width constant and decrease the maximum half-width, the 
sojourn times the particles spends in the bulges decreases causing the mobility to increase 



and approach the maximum value for a straight channel (oj^ 



and e = 1), cf. Fig. \5\ 
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FIG. 5: (Color online) The nonlinear mobility (a), (b) and scaled nonlinear mobility (c),(d) are 
depicted for different scalings of the geometry: constant-bottleneck-scaling, i.e. uj m i n = const, in 
(a) and (c) ; constant-ratio-scaling, i.e. e = w m i n /u; max = const, in (b) and (d). 

(a). In contrast, within the constant-ratio-scaling, where the bottleneck half- width scales 
with the maximum half-width, the nonlinear mobility p depends for small radii only slightly 
on w max , cf. Fig. 0(b). However, with increasing particle size, i.e. radii ratio p, the nonlinear 
mobility p shows a striking dependence on the maximum half-width. 

As pointed out already, with increasing particle radius the effective maximum half-width 
decreases. The effective maximum half-width will exhibit a linear dependence on the radius 
if the particles curvature is larger than that of the channel's boundary function (in physical 
units), 



cff 



-R, for 1/R> 



-d 2 u(x max )/dx 2 



^2^3/2 ' 



(24) 



[1 + (dw(x max )/dx) ]" 

where x max denotes the x-values for which the boundary function assumes a maximum. For 
larger particle radii the effective boundary function shows a kink and the effective maximum 
half-width u)^ ax decreases faster than linearly with the radii ratio p. As the sojourn times 
the particle spend in the channel's bulges depends mainly on the effective maximum width, 
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a nonlinear dependence of the nonlinear mobility is observed for larger particle radii causing 
a peak in the scaled nonlinear mobility, cf. Fig. [5] (c) and (d). 

V. CONCLUSIONS AND OUTLOOK 

We studied the transport of finite Brownian particles through channels with periodically 
varying width. For point size particles it was shown previously js], 0], that the transport 
through such channels could be approximately described by means of the so-called Fick- 
Jacobs equation which is based on the assumption of a fast equilibration in orthogonal 
transport direction. Validity criteria for the capability of this approximation include a 
dependence on the channel shape and predict an upper limit for the force value. In case 
of spherical, finite size particles the maximum force value up to which the Fick- Jacobs 
equation could be applied depends also on the size of the particle. By comparison of the 
approximative result for the nonlinear mobility and the numerical ones we have shown, that 
the equilibration assumptions holds for a wider force-range in case of larger particles than 
it is the case for smaller ones. 

In addition, we pointed out, that the transport of finite, spherical Brownian particles 
in channel geometries with highly corrugated channel walls exhibits some striking features 
which may allow for the development of newly separation devices which extends the func- 
tionality of the sieves. In particular we found, that the nonlinear mobility of Brownian 
particles in such channel structures deviates from the one-over-size dependence predicted 
by the Stokes law for Brownian particles moving in an environment without geometrical 
constrictions. Instead, there is an optimal particle size for which the nonlinear mobility as 
compared to the free mobility exhibits a maximum value. 

Our present study also implicitly used a small concentration of spherical particles such 
that both, effects of particle-particle interactions and forces between particle-particle and 
particle-walls due to hydrodynamic interactions can safely be ignored. These complications 
would require totally new and extensive studies that are beyond this present study. More- 
over, as emphasized with the abstract already, we assumed throughout perfect spherical 
symmetry. Deviations from such spherical symmetry would also impact the viscous friction 



law behavior 



32] and, as well, may give rise to additional, new entropic effects. All such 



complications are beyond the work presented here; all these latter complications, however, 
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open up avenues for interesting future investigations. 
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